Controlled decoherence in a quantum Levy kicked rotator 
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We develop a theory describing the dynamics of quantum Icicked rotators (modelling cold atoms 
in a pulsed optical field) which are subjected to combined amplitude and timing noise generated 
by a renewal process (acting as an engineered reservoir) . For waiting-time distributions of variable 
exponent (Levy noise), we demonstrate the existence of a regime of nonexponential loss of phase 
coherence. In this regime, the momentum dynamics is subdiffusive, which also manifests itself in 
a non-Gaussian limiting distribution and a fractional power-law decay of the inverse participation 
ratio. The purity initially decays with a stretched exponential which is followed by two regimes 
of power-law decay with different exponents. The averaged logarithm of the fidelity probes the 
sprinkling distribution of the renewal process. These analytical results are confirmed by numerical 
computations on quantum kicked rotators subjected to noise events generated by a Yule-Simon 
distribution. 
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I. INTRODUCTION 

The coupling of an open quantum system to an ex- 
ternal reservoir induces a dynamical loss of phase coher- 
ence that results in an irreversible destruction of inter- 
ference effects between states of the system [l|, [2| • This 
environment-induced decoherence efficiently suppresses 
macroscopic superpositions of states and therefore plays 
a central role in our understanding of the transition from 
quantum to classical physics. An important and often 
annoying feature of decoherence for quantum-mechanical 
applications is that it typically occurs very rapidly with 
an exponential decoherence factor, I'(t) = exp(— t/^c), 
where tc is the short coherence time [1, 2]. In some cases, 
however, the coherence time can be controlled by prop- 
erly engineering the reservoir and/or the interaction be- 
tween the system and the reservoir [3] . The feasibility of 
such schemes has been experimentally demonstrated in a 
variety of systems that include cavity QED Q , ion traps 
and cold atom experiments [1, 0, 111 ■ Besides their ob- 
vious technological implications, these experiments pro- 
vide essential tests of the general concept of environment- 
induced decoherence. 

are designed to 
Til , a paradigm 



The atom-optical experiments 
realize the quantum-kicked rotator 
of complex chaotic dynamics which can be simply seen 
as a particle moving on a ring and periodically kicked 
in time. For sufficiently strong kicking strength, the 
classical motion of a kicked rotator is chaotic and the 
evolution in momentum space is diffusive. In the quan- 
tum regime, however, momentum diffusion is asymptot- 
ically suppressed due to dynamical locahzation [12|. [isj. 
In the recent atom-optical reservoir-engineering experi- 
ments, the loss of coherence has been induced by am- 
plitude noise and timing noise [1]. In the first case, 
the amplitude of the periodic kicks is randomly modu- 




lated in a controlled manner, whereas in the second case 
the period between successive kicks is randomly varied. 
In both situations, dynamical localization is extenuated 
leading to a non-vanishing quantum momentum diffu- 
sion with a renormalized diffusion coefficient, in agree- 
ment with theoretical predictions 0, Hi ■ These ex- 
periments enjoy a high degree of control and tunability, 
which makes the atom-optical setup ideally suited to in- 
vestigate the decoherence-induced quantum-to-classical 
crossover of complex quantum systems [l6j . 

One point is especially noteworthy: reservoir engineer- 
ing has so far been used to tune the coherence time, but 
not the exponential time dependence of the decoherence 
factor itself. Our aim in this paper is to describe a 
reservoir coupling scheme for the cold atom experiments 
which allows to modify the time dependence of I'(t) and 
slow down the loss of phase coherence in a controlled way. 
To this end, we superpose the regular periodic pulses with 
noise events occurring with randomly modulated time in- 
tervals r (z. e., we combine amplitude with timing noise). 
We generate these random time intervals by a renewal 
process with a waiting-time distribution w{t). Of spe- 
cial interest is the case where this distribution asymptot- 
ically behaves as a power law, w{t) cx t~^~" . By tuning 
the exponent a of this so-called Levy noise [13], we are 
able to change the mean waiting time between successive 
noise events, given by r = drrw(r), and therefore 
can smoothly interpolate between a fully coupled situa- 
tion where most kicks are perturbed by the noise (large 
a) to an almost isolated situation where most kicks are 
not perturbed by the noise (small a). 

Divergent moments are a hallmark of Levy statistics 
and the mean waiting time r becomes infinite when 
a < 1. In this case the noise is nonstationary and noncr- 
godic [H, [l^ . We show that this type of noise modifies 
the decoherence of the atoms in a striking manner: The 
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exponential time dependence of the decoherence factor 
is replaced by a Mittag-Leffler function [S^l, which 
starts out as a stretched exponential and asymptotically 
decays as a power-law. For this functional form, the co- 
herence time is ill-defined. 

In Ref. [2l| we briefly discussed one manifestation of 
the modified decoherence pattern, namely a subdiffusive 
momentum spreading of the atoms. In this work, we 
present a theory that links the subdiffusive behavior to 
the nonstationarity of the noise. This requires us to ex- 
tend the previous theories of Refs. [3, GE] beyond the 
perturbative short-time regime and to account for the ad- 
ditional timing noise from the renewal process. We also 
discuss the limiting distribution function and the inverse 
participation ratio, which provide additional information 
about the dynamics in momentum space. As alternative 
measures of coherence, we then consider the purity and 
the fidelity of atoms subjected to different realizations of 
the noise, which can be probed in an echo experiment 
(23 , [23j . For nonexponential decoherence, we find that 
the purity initially decays with a stretched exponential, 
which is followed by two regimes of power-law decay with 
different exponents. The decay of the averaged logarithm 
of the fidelity is related to the sprinkling distribution of 
the renewal process. Throughout the paper our findings 
are illustrated by comparison to the results from a nu- 
merical implementation of the kicked rotator. 

The structure of this paper is as follows. Section |TT] de- 
scribes the quantum kicked rotator and the environmen- 
tal coupling scheme of combined amplitude and timing 
noise generated by a renewal process. In Sec. IIIII we de- 
rive the nonperturbative expression for the decoherence 
factor and evaluate it for Levy noise. In the case of non- 
stationary noise with a < 1, the decoherence function is 
nonexponential. Section IIVI discusses how the modified 
decoherence affects the momentum spreading, while Sec. 
IVl contains the discussion of the purity and fidelity. The 
results of this paper are summarized in Sec. IVII 



II. THE LEVY KICKED ROTATOR 
A. Model 

The motion of atoms exposed to pulsed standing-wave 
potentials can be mapped onto a quantum kicked rota- 
tor [Tol . [Tl| . i.e., a point particle moving on a ring 
and subjected to periodical kicks. In suitable units, the 
Hamiltonian of the kicked rotator takes the form 

2 ^ 

H=Y+ KnCoseS{t-n + 0+), (1) 

71— — 00 

where the kicking potential depends on the 27r-periodic 
rotation angle 9 and the kicking amplitude Kn can in 
general be kick dependent. 

From kick to kick, the stroboscopic classical dynamics 



of the kicked rotator is generated by the map 

p{t + l)=p{t) + Kt sin 9, (2) 
9{t+l) = 9{t)+p{t + l). (3) 

Starting from a given initial state the corresponding 
quantum dynamics tpit+ 1) = F{Kt)'ip{t) is generated by 
the Floquet operator 

F{Kt) = exp{-ihr^f/2) exp{-ih~'^Kt cos 9). (4) 

In both cases, the dynamics consists of a kick in which 
the momentum changes by i^tsin^, followed by a free 
rotation in which the rotation angle 9 increases by p. 

In the absence of noise, the dynamics of the kicked 
rotator is entirely controlled by the constant parameter 
Kn = K. For K ^ 5, the classical dynamics is chaotic 
and the growth of the momentum is on average diffu- 
sive with a variance varp(t) ~ Ddt. The classical diffu- 
sion constant is given by Del — 0]. By contrast, 
in the quantum regime, the kicked rotator exhibits dy- 
namical localization, an interference phenomenon akin to 
Anderson localization in disordered solids, which mani- 
fests itself in an exponentially decaying envelope of the 
quasienergy eigenstates F{k) of the Floquet operator 
[12| . As a result, momentum diffusion is suppressed and 
the variance 

varpoW-^*i*[l-exp(-Vi*)] (5) 

saturates after the quantum break time t* ~ D* /fi^ [Tsj . 
The constant D* is of order of the classical diffusion con- 
stant Del, but is subject to quantum corrections [25j . 

The transition between the quantum and the classical 
behavior of the rotator can be induced by subjecting the 
kicked particle to additional random kicks. In the case of 
amplitude noise , the kicking amplitude is explicitly n- 
dependent, Kn = K + km where the perturbations fc„ are 
random numbers. In the presence of conventional timing 
noise the period between the kicks (here set to unity) 
is slightly modulated. 

In this paper we consider an unconventional combina- 
tion of amplitude and timing noise: kicks always appear 
at periodic instances and the amplitude of some of the 
kicks are perturbed. However, the (integer-valued) time 
intervals between the noise events are generated by a re- 
newal process with waiting-time distribution uj{t) (see 
Fig. [T] for a schematic illustration) . The timing of the 
A/''th noise event {N — 1,2,3,.. .) is hence given by 

AT 

tAr = YTn, (6) 

n=l 

where each waiting time t„ is drawn independently from 
the same probability distribution i-u{t). The strength of 
the individual noise events will be characterized by the 
variance kf^ = k of the perturbed kicks, whereas on 
average ktj^ = 0. 
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FIG. 1: (color online), (a) Schematic illustration of the kick- 
ing sequence in the Levy kicked rotator. The period of the 
kicks is set to unity. The timing of the noise is generated 
in a renewal process, which selects the integer-valued waiting 
times r„ from a distribution cij(r). (b) Waiting-time distribu- 
tions of the Yule-Simon form (|24|l for different values of the 
parameter a. The distributions display power-law tails, spec- 
ified in Eq. (|25p . and hence generate Levy noise of variable 
exponent. The lines are a guide to the eye. 



The usual white amplitude noise is obtained for a 
waiting-time distribution lo[t) = 5r,i, so that every kick 
is perturbed. A more intriguing case is a waiting-time 
distribution which asymptotically behaves as a power law 
(Levy noise), 



'w{t) 



CT 



(7) 



where c is a constant. For a < 1 the mean waiting time 
T diverges. Unlike the familiar white amplitude noise, 
Levy noise [l^ generates a nonstationary process with 
unconventional statistical features that are reviewed in 
the following subsection. 



B. Statistical properties of Levy noise 

The properties of noise generated by a renewal process 
can be conveniently quantified by analyzing the num- 
ber of noise events Af{t',t") within an interval [t",t'). 
The number Af{t, 0) is obtained by inverting the ran- 
dom time t_\f, and hence is known as the inverse random 
time of the process 2?!]. On average, the rate of change 
f{t) = dtAf{t,0) of the inverse random time gives the 
probability that there is a noise event at time t. This 



rate is called the sprinkling distribution and is directly re- 
lated to the waiting-time distribution uj{t) via a Laplace 
transformation. 



Jo l-w{u) 



(8) 



We denote throughout the article Laplace-transformed 
functions by a tilde. 

For small values of u, the power-law waiting-time dis- 
tribution ([7|) delivers 



w{u) 1 — UT 

w{u) ~ 1 — 



CTT 



r(Q; + 1) siuTra 



(a > 1), (9) 
(a < 1). (10) 



For a > 1 , it then follows from Eq. ([5]) that the sprinkling 
distribution takes a constant value at large times. 



(a>l). 



(11) 



As a result, the noise asymptotically becomes stationary 
in this case. However, for power-law waiting-time distri- 
butions with a < 1, the sprinkling distribution asymp- 
totically decays over time according to 

f{t) - r~^asin(7ra)/(7rc) (a < 1). (12) 

Due to the exphcit time dependence of f{t), the noise 
here remains nonstationary. We note that nonstation- 
ary also implies nonergodicity, since time averages cannot 
be replaced by ensemble averages [l^. A direct conse- 
quence of Eq. (|12p is that the mean inverse random time 
increases sublinearly. 



J\f{t, 0) ~ sin(7ra)/(7rc) (a < 1). (13) 

More detailed information about the statistical prop- 
erties of the renewal process can be obtained by analyz- 
ing the complete distribution function V{N\ t', t") of the 
number of noise events. The latter is most conveniently 
characterized by the moment-generating function 

oc 

M{z-t',t") = exp[zAA(t',t")] - ^e"^7'(AA;t',t"). 

ji=0 

(14) 

For general time arguments t' and t", this function can 
be expressed as 



M{z-t',t") = M{z-t',0) 



t" 



-(e^-1)/ dsf{s)M{z-t' -s,Q), 
Jo 

(15) 

so that it suffices to study the moment-generating func- 
tion of the inverse random time Af{t, 0). 

We first relate the distribution function V{N;t,Q) of 
the inverse random time to the probability distribution 
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P{tj\f;Af) of the Mth random time, which succeeds via 
the intuitive expression 

V{Af;t,0)= f dt'[P{t';Af)-P{t';Af+l)]. (16) 
Jo 

Since Eq. ^ involves a sum of A/" independent waiting 
times, the Laplace transform of P{tj^;J\f) factorizes, 

P{u;Af) ^ [w{u)]^, 



(17) 



and Eq. p6|) thus yields 



V{Af;u,0) 



1 



dte-''*riz;t,0) = - [l-w{u)][w{u)] 



(18) 

(As indicated, the Laplace transformation is taken with 
respect to the later time argument). By using the defini- 
tion p4| of the moment-generating function Ai{z;t,0), 
we find that its Laplace-transform can be expressed as a 
function of the sprinkling distribution, 



M(z:u,0) = = — . 



(19) 



One should note the explicit periodicity, A4{z;u,0) = 
A4{z + 27ri; u, 0), which is a consequence of the discrete- 
ness of the inverse random time. 

When next going back to the time domain, we again 
have to distinguish between asymptotically stationary 
and nonstationary noise. In the former case (a > 1), 
the moment-generating function is asymptotically expo- 
nential, 

Al(z;i,0) - exp[(e^ - l)t/T] (a > 1). (20) 

Equation characterizes the counting statistics of a 
Poisson process with mean rate t~^. In the nonstation- 
ary case a < I, inserting the asymptotic behavior (jlOp 
produces a moment-generating function of the form 

M{z;t,0) - Ea[t°'{e'--l)T{l+a){simra)/{TTc)] {a < 1), 

(21) 

where Ea is the Mittag-LefHer function |23], defined by 



Eciz) = J2 



-r(an+l) 



(22) 



71=0 



Expression (|2ip becomes exact for all times when one 
passes to continuous inverse times (e^ — 1 —^ z) and con- 
siders the particular case of a scale- free Levy process [l^l ■ 
The moment-generating function then starts out as a 
stretched exponential, M{z;t,0) ~ exp{t"/[r(— a)cic]}, 
while for large times it crosses over to the power law 
A4{z;t,0) ~ {ctc/a)t~°'. For more general waiting-time 
distributions, a reasonable approximation is to adapt the 
time argument by using the mean inverse random time, 
such that 



Miz;t,0)KEa {e' ~l)T{l + a)Af{t,0) (a < 1). 

(23) 

This results in the same power-law asymptotics for large 
times, and accounts for transient behavior in the initial 
decay. 



C. Numerical implementation 

Throughout this paper we compare analytical results 
to results obtained from a numerical implementation of 
the Levy kicked rotator. In the present section, we de- 
scribe the parameters used for these computations. 

Dynamical localization is expected to be the strongest 
when h/{2iT) is approximating a quadratic irrational co = 
In order to render the Hilbert space finite 
we choose a = 24, set Ti = 2ttM/N, and use the third 
convergent of the continued-fraction representation of the 
quadratic irrational, 2l+2l+l4:' This gives M — 577, 
N = 13872. 

The integer N is the Hilbert space dimension, which 
now is large but finite, corresponding to a quantized mo- 
mentum p = hi, I = 0, 1, 2, ... A'^ — 1. Similarly, in po- 
sition space the momentum quantization carries over to 
discretized positions 9 = 2'kI/N. The integer M deter- 
mines the quantum-mechanical periodicity in momentum 
space {p and p + 2itM are equivalent). Each classical 
period of 27r covers N/M « 24 discretized momentum 
values. 

In all our computations, the particle is initially pre- 
pared in the zero-momentum state "0(0) = \p = 0). The 
subsequent propagation is obtained by successive appli- 
cation of the Floquet operator ^ , which can be broken 
down into the application of diagonal matrix multiplica- 
tions in momentum and position representation, inter- 
vened by Fast Fourier Transformations for the passage 
between both representations. 

In each propagation step, the kicking strength is deter- 
mined according to the following procedure. The regular 
kicking strength is set to K ~ 7.5. The noisy perturba- 
tions kn are taken from a uniform box distribution over 
an interval {—W,W), so that k = This noise 

is only applied at times selected by a renewal process, 
which we generate by the so-called Yule-Simon distribu- 
tions (23 



c.(r) 



(24) 



r(r-|-a-|- 1 

These distributions are paradigms of integer- valued prob- 
ability distributions with a power-law tail, 

uj{t) ^ ar{a + l)T-^-°'. (25) 

The mean waiting time is given by 

a 



a-1 



(«>1) 



(26) 



and diverges for a < 1. In the latter case, the asymptotic 
time dependence of the sprinkling distribution is 



7rr(a-|-l) 



(27) 



These properties directly follow from the Laplace 
transform of the waiting-time distribution, which reads 



'w{u) 



a + 1 



e-"2Fi(l,l,a + 2;e-"), (28) 
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FIG. 2: (color online). Time dependence of the variance of 
momentum varpo(i) for a noiseless kicked rotator with regular 
kicking strength K — 7.5. The dashed curve is the theoretical 
prediction ([5]) of dynamical-localization theory [ij . A single- 
parameter fit delivers the value D* = h?t* = 45.28, which is 
used throughout the rest of the present work. The inset zooms 
in onto the region where dynamical localization is established. 



where 2F1 is the hypergeometric function. 

The variance of the momentum for the kicked ro- 
tator in the absence of noise is shown in Fig. [21 A 
single-parameter fit to the prediction ([S]) of dynamical- 
locahzation theory dehvers the value D* = Ti^t* = 45.28, 
which is our only fitting parameter and will not be ad- 
justed any further throughout the rest of this work. 



III. DECOHERENCE IN QUASIENERGY 
SPACE 



The effect of noise on the quantum kicked rotator has 
been first investigated by Ott, Antonsen and Hanson, 
who developed the following intuitive picture (14i | . In the 
absence of noise the quantum system is fully coherent: 
the quasienergy eigenstates of the Floquet operator F{K) 
are exponentially localized in momentum space, chaotic 
diffusion is strongly suppressed and the momentum diffu- 
sion constant asymptotically vanishes. The effect of the 
external noise is to couple the quasienergy states and to 
induce transitions between them. As a result, dynamical 
localization is extenuated and quantum diffusion takes 
place. 

A systematic approach for calculating the non-zero 
quantum diffusion constant in the presence of station- 
ary noise has been developed by Cohen 'l5i|. It is based 
on an analysis of the decay of a decoherence factor which 
is related to the survival probability of the quasienergy 
states of the Floquet operator. In Ref. 15| this deco- 



In Ref. 

herence factor was calculated perturbatively for short 
times (the emphasis was on stationary noise with ar- 
bitrary correlations, including the intricate correlations 
of zero-temperature noise). In the following, we use 



a simple random-phase approximation to calculate the 
decoherence factor beyond the perturbative short-time 
regime. This approximation is valid for random uncorre- 
lated noise events and remains immediately applicable to 
noise generated by a renewal process, including the case 
of nonergodic, nonstationary Levy noise. 



A. Quasienergy transitions and decoherence factor 

The quasienergy eigenstates \r) of the noiseless system 
are defined by the eigenvalue equation 



F{K)\r) =exp(-ie,)|r). 



(29) 



where is the quasienergy. A noise-free propagation 
step does not alter the quasienergy state, while a noisy 
kick introduces transitions with amplitude 

{s\F{Kn)\s') — exp(— ies)(s| exp(— i?i^"^A:„ cos0)|s'), 

(30) 

where s' and s are the quasienergy index before and after 
the propagation step, respectively. 

Over a larger time interval [i",i'), the survival ampli- 
tude in a quasienergy state \r) is defined as 



Ar{t'X) = (r|r [] F(X„)|r)exp[z(i'-t")e.], (31) 



where T is the time-ordering operator. In this expres- 
sion, we have explicitly provided for two time argu- 
ments, as is required when working with a renewal pro- 
cess which may produce nonstationary noise. The ex- 
ponential factor compensates for the dynamical accumu- 
lation of quasienergy phases in the noiseless system, for 
which Arit' ,t") = 1. 

The averaged survival amplitude is denoted by 



(32) 



where the average is over the index r as well as over the 
timing and amplitude of the noise. 

The decoherence factor Vit' ^t"^ that later appears in 
the momentum spreading and in the decay of the purity 
is given by 



V(t',t") = A,(t',t")Al{t',i"), 



(33) 



where the average now extends over both quasienergy 
indices r ^ s. 



B. Average over amplitude noise 

In order to analyze the decoherence factor, we first as- 
sume that the timing of noise events is fixed and average 
over the random detunings fc„ of the noise events, as well 
as over the quasienergy indices. These averages become 
straightforward within a random-phase approximation. 
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which uses that the random detunings are uncorrelated 
and moreover exploits features of the complex quantum 
dynamics encoded in the quasienergy transition ampli- 
tudes ([50]) . For this we have to distinguish amplitudes 
which preserve the quasienergy index (diagonal matrix 
elements, which are of modulus one if there is no noise 
and less than one in the presence of noise) from ampli- 
tudes describing transitions between quasienergy states 
(off-diagonal matrix elements, which are only finite in the 
presence of noise) . 

For noise strength k, the averaged modulus of the 

diagonal matrix elements in Eq. (j30p can be obtained in 
a simple expansion, 



exp{-in~^ki cos 61) w 1 - {K/2h^) cos2 9 = 1- k/AH^, 

(34) 

which delivers the survival time 



(35) 



We will shortly see that tc coincides with the decoher- 
ence time for conventional amplitude noise, which acts 
stationarily in every kick. The condition set for the noise 
strength implies tc ^ 1- This assumption is not very 
restrictive, since it only prevents total decoherence by a 
single event, but does not impose any restriction on tc 
with respect to the localization time t* ~ K'^/2h'^, which 
is also much larger than unity. We hence reserve the 
notion of weak noise for the more restrictive condition 
K <C h'^/K'^, which implies tc'^ t* . 

Each noise event mixes the initial quasienergy state 
into about (, = t* 1 states whose localization centers 
are in reach of the localization length [25| . Each individ- 
ual off-diagonal transition amplitude hence is small, and 
for weak noise {tc >• t*) the averaged off-diagonal matrix 
elements indeed can be trivially neglected. The averaged 
Floquet operator F{Kn) is then simply approximated by 



(s|F(ir„)|s') -5..s'diag[e- 



-i/(2te)i 



(36) 



which is diagonal in quasienergy space. 

When tc and t* are of comparable order, individual 
contributions which include quasienergy transitions arc 
still small, but in the composed survival amplitude (|3ip 
one cannot disregard the proliferation of such terms when 
the number of transitions is increased. It then becomes 
necessary to exploit the consequences of the complex 
dynamics of the kicked rotator, which induces random 
phases into the off-diagonal transition matrix elements. 
The random phase has several separate origins: 

(1) For fixed s and s', the transition amplitude ((30|) 
sensitively depends on the random detuning fc„ of the 
kick. 

(2) In momentum space, the quasienergy states have 
an exponentially decaying envelope, but a quasi-random 
internal structure. Consequently, even after the average 
over the detuning, the phase of the transition amplitude 
depends strongly on the quasienergy indices s and s' . 

In the ensuing random-phase approximation, the phase 
of the off-diagonal transition amplitudes is assumed to be 



uniformly distributed in [0, 27r), while the squared mod- 
ulus of these matrix elements is of order 2/(tct*) (corre- 
sponding to the bandwidth f = i*/2). The squared mod- 
ulus of the diagonal elements is given by exp(— 1/ic) ~ 

1 - l/tc- 

In this approximation, the averaged Floquet operator 
F{Kn) simply retains its diagonal form ([55)1 . This is 
all we need for the computation of the averaged survival 
amplitude V^-^^t' ,t") and leads to the intermediate ex- 
pression 

V^'^ (f, t") = eM'Mit',t")/{2tc)], (37) 

where Af(t' , t") is the inverse random time of the renewal 
process as introduced in Sec. HH 

The random phase approximation also establishes a di- 
rect connection between V'-^^t' ,t") and the decoherence 
factor I?(t', t"). Since in Eq. ([55)) the two quasienergy in- 
dices r and s are not identical, there are no contributions 
which link the quasienergy states in exactly the same se- 
quence. The decoherence factor therefore immediately 
factorizes into the averaged survival amplitudes of the 
two quasienergy states, yielding 



V{t',t") = [I?(i)(t',i")]2 = exp[~J\f{t',t")/tc]. (38) 

The remaining average over the timing of the noise events 
does not factorize, but because of the simple exponential 
form of Eq. ([57)) . merely amounts to a factor of two in 
the exponent. We now analyze this average for different 
types of renewal processes. 



C. Average over timing noise 

In the conventional case of stationary noise with 
lij{t) = (5t-,1i where timing noise is absent, the inverse 
random time, Af{t' , t") = t" — t', does not fluctuate. Ac- 
cording to Eq. ([38)) . the decoherence factor 



V{t', t") = exp[-(t' - t")/tc] (stationary noise) (39) 

then decays exponentially over time, and tc is identical 
to the coherence time. 

For noise generated by a renewal process with power- 
law waiting-time distributions, the subsequent analysis 
depends on whether the process is characterized by a 
mean waiting time r, or shows the nonstationary fea- 
tures of Levy noise. This bifurcation can be made explicit 
by realizing that Eq. (|37p is identical to the moment- 
generating function A4{z;t' ,t") of the number of noise 
events, introduced in Eq. p^ . which has to be evaluated 
at the noise-strength dependent value z = —l/tc- We 
therefore can immediately exploit the results of Sec. Ill Bl 
(Since we assume tc ^ 1 we can set — 1 « —l/tc ) 

It follows from Eq. that, when the noise is asymp- 
totically stationary {a > 1), the decoherence factor re- 
tains its exponential form. 



V{t',t") - exp[-(t' - t")/Ttc] {a > 1). 



(40) 
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FIG. 3: (color online). Averaged distribution of momenta P{p;t) at selected times t for kicked rotators which are subjected 
to Levy noise generated by a Yule-Simon distribution. The average is over 10^ realizations. The smooth curves are fits to 
double-sided exponential and Gaussian distributions. The unit of momentum is the localization length p* — D* /Ti. (a) In the 
upper panels a = 2.0, k = 1/30000, corresponding to a decoherence time tc — 12.3 1* (where t* is the localization time). At 
i = 15t* « 1.2 tc the distribution function fluctuates around a double-sided exponential function. At t = 30t* ~ 2.4 the 
distribution function is still exponential, but the fluctuations are suppressed. At t — 120 1* « 9.8 the distribution function 
is well described by a Gaussian, (b) In the lower panels a — 0.5, n — 1/3, corresponding to stronger but nonstationary, 
increasingly rare noise events. In this case, a decoherence time cannot be defined and the distributions function remains 
exponential for all times t > t* . Note that the peak at p = arises from coherent backscattering. 



The effective coherence time Ttc is directly proportional 
to the mean waiting time of the noisy kicks and therefore 
increases with decreasing exponent a. In particular, the 
effective coherence time becomes infinitely large when a 
drops below unity. According to Eq. ((23)) . the functional 
form of the decoherence factor then changes to 



V{t,0)^ Eo, -T{l + a)J\f{t,0)/t^ (a<l). (41) 



After a (typically short) transient, J\f{t,0) ~ 
sin(7ra)/(7rc). Hence, we have the remarkable result 
that Levy noise induces a nonexponential loss of coher- 
ence which typically starts out as a stretched exponential, 

X'(t,0) ~ exp[r/(r(-a)cic)], (initial decay), (42) 

and for large times crosses over to a power law, 

I?(t, 0) ~ (ctc/a)i~" (asymptotic decay). (43) 

The functional dependence of the decoherence factor can 
thus be changed in a controlled manner by tuning the 
value of the exponent a. We now explore how this affects 
the observable properties of the kicked rotator. 



IV. MOMENTUM SPREADING 

In the atom-optical experiments, the momentum can 
be observed by a time-of- flight measurement, where the 
kicking sequence is terminated at a specified time t and 
the atom cloud is allowed to expand ballistically to a 
much larger size. A snapshot is then taken and the travel 
distance is converted into the momentum. This technique 
provides direct access to the complete momentum distri- 
bution P{p; t). 

In the presence of dynamical localization, the mo- 
mentum distribution shows the exponential envelope of 
the atomic wavefunctions in momentum space, while a 
Gaussian envelope is expected when diffusive momen- 
tum spreading is induced by stationary timing or am- 
plitude noise. In order to give an early indication of 
the effects of nonexponential decoherence, we compare 
in Fig. [3] the momentum distribution of Levy kicked ro- 
tators with asymptotically stationary and nonstationary 
noise {a — 2 and a — 1/2, respectively). For a — 2, 
the strength of the noise is set to k = 1/30000, corre- 
sponding to a decoherence time tc — 12.3 i*. For t ^ tc 
the distribution function is well described by a double- 
sided exponential function; when the time is increased 
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FIG. 4: (color online). Same as the left panels in Fig. O but 
for K = 1/300 [tc = 0.12f*), so that t = 15t* > U. 

it eventually evolves into a Gaussian. This follows the 
expectations for conventional stationary noise. By con- 
trast, for a — 1/2, where a well defined coherence time 
does not exist, the probability density keeps its double- 
sided exponential shape at all times. This is the case 
even though we increased the strength of the individual 
noise events by four orders of magnitude to k = 1/3, so 
that tc — 0.00123 is much reduced. This indicates that 
diffusive momentum spreading is never attained, even in 
the limit of very strong noise events. 

Figure m shows the momentum distribution dXt = IM* 
for the more moderate noise strength k, = 1/300 {tc — 
0.12 1*) which is used in several subsequent figures in this 
work. For a = 2.0, the momentum distribution is always 
well-approximated by a Gaussian, while for a = 0.5 and 
t > i*, it is well-approximated by a double-sided expo- 
nential function as soon as t > t* . 

We now proceed to characterize the evolution of the 
momentum distribution by two quantities, the variance 



varp and the inverse participation ratio (IPR). 

A. Variance of momentum 

The temporal evolution of the momentum distribu- 
tion P{p] t) has two origins. The first component is the 
spreading of the momentum wavefunction of an individ- 
ual atom with respect to the mean momentum of the 
atom. This momentum spreading is characterized by the 
variance. 



varp(t) = (p2) - (p)2. (44) 

As indicated, this variance is computed from the 
quantum- mechanical expectation values, and then aver- 
aged over the noise. The second component of the evo- 
lution of the momentum distribution is the spreading of 
the mean momenta of the different atoms. Hence, in 
general, the variance is smaller than the total vari- 
ance computed from the momentum distribution P{p; t). 
In the particular case of the kicked rotator, however, 
this difference becomes negligible, due to the symmetry 
(0,p) {—9,—p) which holds even in presence of the 
noise. 

In the following, we first establish a general relation 
between the variance of the momentum of the noise-free 
and noisy systems, and then evaluate this relation for 
the Levy kicked rotator. We will then find that nonex- 
ponential decoherence manifests itself in a subdiffusive 
spreading of the momentum. 

As in the earlier theories of delocalization due to de- 
coherence [15], the starting point of our considerations is 
the force-force correlation function. 



C(t',i") = {Kt'Kr, sm9t> sin^t.), (45) 

which is composed of the momentum increments in an 
individual time step, as given by Eq. ([2]). In the above 
equation, 9t is the angle operator in the Heisenberg pic- 
ture. The average includes an average over the initial 
momentum of the particle, so that the expectation value 
is replaced by the normalized trace A'^~^tr ^v, where N is 
the Hilbert space dimension. 

The total change of the momentum is obtained by inte- 
grating the force over time. The variance therefore reads 

t-i 

varp(t)= C{t',t"). (46) 

t',t"=0 

In order to establish the relation to the decoherence 
factor, the force-force correlation function is expanded in 
the basis of quasienergy eigenstates [11] , 
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ri+l)*{Sm+l\F{Km)\Sm,)- (47) 

{rn,Sn} l.m—t" 



The prime in the sum over quasienergy indices excludes 
the initial index rt" , which originates from the normal- 
ized trace and hence instead is averaged over. 

We first apply the random-phase approximation to the 
matrix elements of the sine kicking potential, which in 
principle selects three types of terms: (i) terms with = 
Tt", Sf = St", but Tf ^ Sf] (ii) terms with rf = s^, 
ft" = St" , but Tf 7^ Tt" ; (iii) terms with = rt" — st' = 

St". 

The terms of type (ii) and (iii) involve the expec- 
tation value of the force, which vanishes in the spe- 
cific case of the kicked rotator, thanks to the symmetry 
{0,p) {—d,—p) (a careful analysis of such terms will 
be required for the purity, see Section |V|. For the terms 
of type (i), on the other hand, the product of transition 
matrix elements {ri\F(Ki)\ri-j-i) in quasienergy space ex- 
actly arranges itself into the definition of the decoher- 
ence factor This correspondence also holds when 

C{t',t") is evaluated for t' < t" , provided we define 
V{t',t") =V{t",t'). 

For t' = t" , the average over the noise also gives rise to 
a contribution to the typical force acting on the particle 
of the form 'K^{sm9^) = + {K/2)f{t'), where f{t) 

is the sprinkling distribution. Collecting all terms, we 
arrive at 

C{t',t") = Co{t',t")V{t',t") + ^f{t')St'.t", (48) 

where Ca{t',t") is the force-force correlation function in 
the absence of noise. Besides the correction from the 
sprinkling distribution, which is small for weak noise, this 
recovers the result from the perturbative short-time the- 
ory in Ref. ■ It follows that the correct extrapolation 
of the decoherence factor beyond the initial decay is given 
by Eq. 

The variance of the momentum follows when Eq. ([^5]) 
is inserted into Eq. (|46p . It is then convenient to pass 
to a continuous time argument and perform integrations 
by parts, so that all force correlations are expressed in 
terms of the variance of the momentum in absence of the 
noise, which we denote by varpo(0- then obtain as 
one of our main results the following expression, 

varp(t) = vSirpo{t)V{t,0) + ^7V(t,0) 

+ / ds var pQ{t — s) ds'D{t, s) 
Jo 

— / dsvarpols) dsV{s,0) 
Jo 

-f ds' f ds"v(iipo{s' - s")ds'ds"V{s',s").{49) 
Jo Jo 



I 

The second term on the right hand side can be neglected 
for small noise strength k. Equation (j49p expresses the 
momentum variance in presence of noise exclusively in 
terms of the variance in absence of the noise, and deriva- 
tives of the decoherence factor T?{t',t"), which play the 
role of memory kernels. Because of the monotonicity 
properties of 'D{t',t") all terms are positive. Since the 
result is expressed directly in terms of varpo, Eq- (011) 
moreover accounts for all dynamical correlations of the 
noiseless dynamics. In particular, the random-phase ap- 
proximation used in the derivation only relies on the com- 
plex quantum-dynamics of the noiseless system, but does 
not require that it displays dynamical localization. 

In exploring the consequences of Eq. (|49p for the kicked 
rotator, we combine Eq. ([5]) for the momentum variance 
in absence of the noise with the various decoherence fac- 
tors found in Section [1111 First, let us again make contact 
with the known results for stationary noise, for which the 
decoherence factor decays exponentially as given in Eq. 
pop. The corresponding momentmn spreading is 

""'^^ = TTTjF^ + (1 + J/t,r "P^"'''** 

_ (50) 

This result was first proposed in Ref. I15| as one of two 
possible extrapolations of the perturbative short-time re- 
sults, and has been successfully applied to quantify the 
decoherence process in the experiment of Ref. Q. The 
present derivation shows that the other proposed extrap- 
olation indeed was rightfully discarded. For large times, 
the second term approaches a constant and can be ne- 
glected in comparison to the first term, which increases 
linearly in time. Hence, for stationary noise when quan- 
tum coherence is lost exponentially, the momentum dif- 
fuses asymptotically with a renormalized diffusion con- 
stant D*/{l + tc/t*). 

In the presence of asymptotically stationary noise gen- 
erated by a renewal process with a > 1, the decoherence 
function is given by Eq. (I40p . in which the coherence 
time is modified by the mean waiting t. This does not 
induce any qualitative change in the momentum spread- 
ing, which is still diffusive according to a renormalized 
diffusion constant D* /{I + Ttc/t*). 

As a approaches unity from above, the decoherence 
time increases, and the quantum diffusion constant van- 
ishes at a = 1. It is clear that this has to be accompanied 
by a qualitative change of the momentum spreading it- 
self. For a < 1, the asymptotic spreading is dominated 
by the last term in Eq. (|49p . which can be approximated 
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varp(i) ~ — varpo(oo) / ds' i ds" ds'ds"'D{s' , s") 
Jo Jq 

r)*f* 

= ^U{t,0). (51) 

Using the power-law asymptotics (jl3p of the inverse ran- 
dom time for a < 1, we immediately find that the mo- 
mentum spreads subdiffusively, 



, . D*t* sinTra „ 
varp(i) 1" 

tr TTC 



(52) 



The full time dependence of the momentum variance 
of Levy kicked rotators with asymptotically stationary 
or non-stationary noise is shown in Fig. [5l The upper 
panel shows the results for a = 2, for which the momen- 
tum spreading is diffusive. The classical diffusion is at- 
tained for moderate noise strength k > 1/300. The lower 
panel shows the results for a = 1/2, and confirms that 
in this case the spreading is subdiffusive, varp cx t-^^^. 
Note that this behavior persists even when the strength 
of the individual noise events is large. The theoreti- 
cal predictions in the plots follow from Eq. (|49p . where 
varpo(i) is taken from Eq. (O, while 'D{t',t") is deter- 
mined for Eq. psp [i.e., the moment-generating function 
of J\f{t',t") for the Yule-Simon distribution is evaluated 
at z = —1/tc — ~Kj2Ti'\. There is good agreement be- 
tween the numerical results and the theoretical predic- 
tions. 
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FIG. 5: (color online). Time dependence of the momentum 
spreading varp(f) for kicked rotators which are subjected to 
Levy noise generated by a Yule-Simon distribution with noise 
exponent a = 2.0 (panel a) and a = 0.5 (panel b). The results 
of numerical computations (solid curves) are compared to the 
theoretical prediction from Eq. (|49p (dashed curves). The 
straight solid line labelled "classical" is the classical diffusion. 



B. Inverse participation ratio 

In this section we study an alternative measure of 
the momentum spreading, the inverse participation ratio 
(IPR), and show that this quantity also reveals the exis- 
tence of two distinct (exponential and non-exponential) 
decoherence regimes. The inverse participation ratio is 
defined as 

IPR(t) ^^hMiJF, (53) 

n 

where '0n are the amplitudes of the wavcfunction in a 
fixed basis. The IPR can be interpreted as a measure for 
the inverse dimension of an effective Hilbert space ex- 
plored by the quantum dynamics. In general, the IPR 
is basis dependent. In the following, we provide a quasi- 
classical estimate of the IPR by introducing basis states 
which cover Planck cells of area 27r?i in phase-space and 
are centered at phase-space positions (0„,p„). This esti- 
mate turns out to be independent on the details of the 
quasiclassical basis. 

For a given time, the complex amplitudes fluctuate 
as n is varied, corresponding to speckle noise in the phase- 
space position {0m Pn)- In order to quantify these fluctu- 
ations, we introduce a phase-space function P{6n,Pn',t) 



such that the statistical average of the squared amplitude 
is 



\^^{t)\^ = 27rhP{0,„p,,;t). (54) 

By assuming that the real and imaginary parts of ipn are 
Gaussian distributed and independent, we have 

W^Fj\^^2WJtW'- (55) 

In the quasiclassical approximation, the function 
P{9,p;t) is now identified with a smooth, properly nor- 
malized phase space density with continuous variables 9 
and p. By combining Eqs. and ([55)) . the IPR can 
then be written in the form 

IPR(t) = / ^ X 2[27rnP{0,p;t)f . (56) 

The above expression is invariant under canonical trans- 
formations of the phase-space variables and hence ex- 
hibits the advertised independence of details of the quasi- 
classical basis. 

In the Levy kicked rotator, the quasiclassical proba- 
bility density becomes ^-independent after a short time. 
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so that to a good approximation P{9,p;t) ~ P{p;t)/2TT. 
This leads to a compact expression for the IPR, 



IPR(<) = '^^ J dpP^{p]t). 



(57) 



The integral in Eq. ([57|l can be calculated in two 
simple cases: For a Gaussian profile, P{p\ t) = 
(27rvarp(i))~^/^ exp[— p^/2varp(t)], typical for diffusive 
spreading, we obtain 



IPR(i) 



n 



TTvarp{t) 



(Gaussian profile). 



(58) 



On the other hand, for a double-sided exponential prob- 
ability distribution, P{p;t) — (A/2) exp(— A|p|) with A — 
•\/2/varp, we find 



IPR(t) 



•\/2varp(i) 



(exponential profile). (59) 



Interestingly, in both Eqs. ([58| and (j59p the IPR is pro- 
portional to ?i/ y^varp{t). In view of Eq. ([52]) . we hence 
expect that the asymptotic decay of the IPR is a power- 
law with exponent —1/2 for a > 1, while for a < 1 
the exponent is reduced to —a/2. These predictions are 
confirmed by the numerical results shown in Fig. [S] We 
can therefore conclude that the quasiclassical expressions 
([55)1 and ([55]) offer a good quantitative estimate of the 
IPR of the Levy kicked rotator. 
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V. PURITY AND FIDELITY 

A convenient quantity to further elucidate the loss of 
coherence in a quantum-dynamical system is the fidelity 

where the states ^p{t) and ip'{t) start out from the same 
initial condition "0(0) — V''(0), but are propagated with 
different realizations of the noise. 

A remarkable property of the fidelit y is that it can be 
directly probed in an echo experiment [3l|. In the latter, 
the quantum state is first let to evolve up to a certain 
time t at which the dynamics is reversed, e.g., by ap- 
plying a pulse that inverts the momentum (j> — > —p). 
When at time 2t the momentum-reverting pulse is ap- 
plied again, the system approximately returns into its 
initial state. The fidelity characterizes the quality of this 
echo, which is never perfect because of imperfections in 
the momentum-reverting pulse and the noise in the dy- 
namics. The destruction of phase coherence in a quan- 
tum system such as the kicked rotator can therefore be 
directly studied in echo-like experiments. 

In the present section, we study various statistical as- 
pects of the fidelity for the case that its decay is induced 
by Levy noise. 



FIG. 6: (color online). Time dependence of the inverse par- 
ticipation ratio IPR for kicked rotators which are subjected 
to Levy noise generated by a Yule-Simon distribution with 
K = 1/300. The main panels are single- logarithmic while 
the insets show the same data in a double-logarithmic pre- 
sentation, (a) Noise with exponent a — 2.0, corresponding 
to a decoherence time tc — 0.12 (where t* is the local- 
ization time) so that the momentum distribution P{p; t) is 
always well-approximated by a Gaussian. The dashed line is 
the quasi-classical prediction ((58}. (b) Noise with exponent 
a = 0.5, where for t > t* the momentum distribution P{p; t) is 
always well-approximated by a double-sided exponential func- 
tion. The dashed line is the quasi-classical prediction (|59p . 

A. Purity 

The statistical average of the fidelity 

Tirj ^trp^it) (61) 

is formally equivalent to the so-called purity of the statis- 
tical mixture of the perturbed states, which is described 
by the density operator 

Pit) = \m){m\- (62) 

The average is over the amplitude and timing noise gen- 
erated by the renewal process. The initial state is fixed. 
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so that the purity starts at tr (0) = 1 . For completely 
incoherent superpositions the purity tends to the value 
1/A^, where is the Hilbert space dimension. 

In the following, we show that the initial and inter- 
mediate decay of the purity is directly related to the de- 



coherence function T>{t). For the long-time asymptotics, 
we moreover establish a connection to the inverse partic- 
ipation ratio. 

In the quasienergy basis, starting from Eq. (|6ip . the 
purity can be expanded as 



t^p\t)= ^r,rr'Jtoi^s'„ n i^w\FiKi)\ri){r[,+,\FiKl,)\r[,)*{sra+^^^^ 

{fn-.Sn.r'^^.s'^} l,m,l' ,m'—0 

(63) 



where the prime in the sum denotes the constraint rt — 
r'^, St = s'^. Furthermore, V'r is the initial wavefunction, 
expanded in the quasienergy eigenbasis, and the values 
Ki and K'^ are obtained from independent realizations of 
the renewal process, which differ both in the timing as 
well as in the detunings of the noisy kicks. 

We first apply the random-phase approximation to the 
expansion coefficients tpr, which selects three groups of 
terms: (i) terms with ro = rg, sq = Sg, but rg 7^ sq; (ii) 
terms with rg = sg, rg = s'q, but rg 7^ rg; (iii) terms with 
ro = rg = So = s'q. 

For each class of terms, the random-phase averages 
over the detunings and timing of the noisy kicks can be 
carried out following the procedure which we applied to 
the decoherence factor and the survival probability in 
quasienergy space, described in Sectionlllll For the terms 
of group (i) , this procedure selects only diagonal matrix 
elements with fixed r„ = r^ = r, s„ = = s, just as en- 
countered for the decoherence factor. Since we now deal 
with two independent renewal processes, they deliver a 
factor exp[—Af{t,Q)/tc + M'{t,0)/tc] which depends on 
the sum of inverse random times in the propagation of 
the two states. After averaging over the timing noise 
of the two independent renewal processes, these terms 
therefore deliver a factor T>^{t, 0). 

Among the terms of group (ii), the random-phase ap- 
proximation requires to pair all amplitudes to probabil- 
ities, r„ — s„, r'^ = s^. In the calculation of the vari- 
ance varp(i), these terms could be neglected because they 
were multiplied by the vanishing expectation value of the 
detuned force. We shall denote these contributions by 

Sr^s,cl{t, 0). 

Among the terms of group (iii), both terms which are 
completely diagonal (type i) as well as terms which are 
paired to probabilities (type ii) survive. A neat sepa- 
ration of these terms follows when we require that for 
r = s, the terms are called of type (ii) when they contain 
at least one quasienergy change. Appropriately attribut- 
ing these terms to the groups (i) and (ii), the purity then 
takes the form 

tr p2 (i) = ^ I n2)2 (^t, 0) + Sr^sAt, 0)] , (64) 



where the restriction r 7^ s is now lifted. 

For the term proportional to I?^(i, 0) we can now ex- 
ploit the normalization of the initial wavefunction. Look- 
ing at the second term, it easy to see that these contri- 
butions can be neglected for short times t <^tc, since per 
definition they contain at least one quasi-energy transi- 
tion. The significance of these terms becomes apparent 
when one re-evaluates the purity for long times. For this 
we start from the general expression, 

tr p2 (<) ^ J2 rnitWnitWmi'^it). (65) 

nm 

and adopt the quasi-classical basis in which we computed 
the IPR (see Section |IVB|) . For large times t ^ tc ''^t 
is reasonable to assume that the relative phases of the 
wavefunction expansion coefficients ipn and ip'^ in this 
basis are completely randomized. Assuming furthermore 
IV'nP ~ IV'nPi we find that trp^(i) asymptotically ap- 
proaches the IPR as defined in Eq. ([55)1 . The terms of 
group (ii) are therefore associated to the eventual satu- 
ration of the purity at the inverse of the effective Hilbert 
space dimension. We hence arrive at our main result for 
the purity, 

tvp^{t)f=iV^{t,{)) + WR{t). (66) 

The initial and intermediate decay of the purity is 
given by the first term, which is simply the square of the 
decoherence function 2?(i, 0). For stationary noise the 
initial decay is thus exponential, tr/5^(t) = exp(— 2t/tc). 
A similar behavior is predicted for Levy noise with a > 1, 
for which tc replaced by rtc- For nonstationary noise 
with a < 1, Eq. ([4T|) predicts that the initial and in- 
termediate decay shows the functional features of the 
squared Mittag-Leffler function, i.e., a stretched expo- 
nential which for later times crosses over to a power law 
with exponent ~-2a. 

The long time behavior of the purity is given by the 
saturation at the inverse effective Hilbert space dimen- 
sion, which we estimate in terms of the inverse participa- 
tion ratio. According to Eqs. and this term is 
proportional to [varp(t)]~^/^. For a > 1 the exponential 




0.001 0.01 0.1 1 10 100 



tif 

FIG. 7: (color online). Time dependence of the purity tr 
of Levy kicked rotators for different noise exponents a and 
fixed noise strength «; = f/300. The results of numerical sim- 
ulations (solid curves) are compared to the theoretical predic- 
tions from Eq. (|66p (dashed curves) 



decay therefore asymptotically crosses over to a power- 
law decay with exponent —1/2. As discussed before, for 
a < 1, 'D^{t) itself already crosses over from a stretched 
exponential to a power law with exponent —2a. As a con- 
sequence of the subdiffusive momentum spreading, this 
intermediate power law is asymptotically replaced by a 
slower power law with exponent —a/2. 

As a function of the inverse random time J^{0,t), the 
crossover to the asymptotic power-law of the saturation 
background sets in at Af{0,t) w tc In t*. For a > 1, this 
condition translates into a well-defined crossover time t w 
Ttc^nt*. For a < 1, however, the nonergodicity of Levy 
noise entails that the crossover is blurred. 

The numerical results presented in Fig. [7] confirm the 
quantitative and qualitative predictions of Eq. (|66|) for 
the purity of the Levy kicked rotator. For a > 1 the 
initial decay is exponential and the asymptotic decay is 
algebraic with exponent —1/2. The background contri- 
bution is here clearly visible and displays the predicted 
power-law decay. For a < 1 the decay is a stretched expo- 
nential which crosses over to an algebraic decay with ex- 
ponent —2a. Because of this significant slow-down of the 
decoherence, the eventual crossover of the exponent to 
—a/2 is only observed in the numerical data for a ^ 0.7; 
for smaller values of a, the background contribution is 
hardly discernible. 



B. Averaged logarithm of the fidelity 

For a system with complex quantum dynamics such 
as the kicked rotator, it is often observed that the fi- 
delity fluctuates significantly from realization to realiza- 
tion [2^ , and from Fig. [7] it is clear that this tendency 
is exacerbated for nonergodic noise. These fluctuations 
are suppressed when one instead considers the averaged 
logarithm lnT{t) of the fidelity. 
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FIG. 8: (color online). Averaged logarithm of the fidelity of 
Levy kicked rotators for different noise exponents a and fixed 
noise strength k = 1/300. The results of numerical simula- 
tions (solid curves) are compared to the theoretical predic- 
tions from Eq. (f55)l (dashed curves), which does not account 
for the saturation background. 

By passing from the fidelity to its logarithm, the multi- 
plicative suppression of coherence per noise event is con- 
verted into additive contributions. For initial and inter- 
mediate times, the averaged logarithm hence probes the 
mean number of noise events in both noise realizations, 



AA(0, t) + AA'(0, t) = 2 / fit') dt', (67) 
Jo 

where f{t) is the sprinkling distribution. This leads to 
the prediction 

2 /■* 

InJF(i) = / f(t')dt' (short and intermediate times), 

tc Jo 

(68) 

while asymptotically a saturation background should set 
in (we do not give an estimate of this background). 

The validity of this prediction is confirmed by the nu- 
merical results in Fig. [8l which shows the averaged loga- 
rithm in a representation that makes it readily compara- 
ble to the results for the purity (Fig. [7]). There is good 
agreement for small and intermediate time, until satura- 
tion sets in, which is not incorporated in Eq. ((S5|). Two 
noteworthy observations in the comparison to the purity 
are first the inequality exp(lnJ^) < trp^, and secondly 
the fact that the statistical fluctuations arc now distinc- 
tively suppressed. This indicates that the logarithm of 
the fidelity probes typical noise realizations, while the 
purity is influenced by rare noise realizations in which 
coherence is exceptionally well preserved. 

VI. DISCUSSION AND CONCLUSIONS 

In this concluding section, we would like to situate our 
work in the general context of the study of anomalous 
diffusion and nonexponential relaxation. The investiga- 
tion of diffusion in complex environments enjoys a long 
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history starting with the seminal work of Scher and Mon- 
troU on dispersive transport in amorphous semiconduc- 
tors [s^ . Subsequently, extensive work on non-Brownian 
motion and nonexponential decay induced by Levy noise 
in space and/or in time in classical systems has been 
carried out, most prominently using the powerful Con- 
tinuous Time Random Walk (CTRW) formalism (for a 
review see Ref. JBSJ). In the last few years, this line of re- 
search has been extended to quantum systems. The novel 
aspect of these developments is the possibility to study 
the interplay of complex dynamics and quantum phenom- 
ena. Thus the interaction with a composite environment 
with extra degrees of freedom, and the entanglement of 
the quantum system with the latter, has for instance be 
shown to lead to power-law decay {34|. The relaxation 
of a quantum two-level system subjected to stationary 
and nonstationary pow er-law noise has been respectively 
examined in Refs. [3^ and [HI. In the latter work, the 
presence of aging dephasing has been demonstrated. Ad- 
ditionally, anomalous fast decoherence induced by spatial 
Levy noise, stemming from a chaotic random-matrix en- 
vironment, has been investigated in Ref. JiT]. 

Our approach in the present work has been to use 
the quantum kicked rotator as a powerful tool to ex- 
plore decoherence and its interplay with complex quan- 
tum dynamics, in particular, nonstationarity. The latter 
property is well-known from the physics of disordered or 
glassy materials [38] . Taking advantage of the unique de- 
gree of tunability of the atom-optical realization of the 
kicked rotator, we have put forward a way to engineer a 
complex reservoir, providing full control over the station- 
arity of the environment or the absence thereof. Specif- 
ically, we have proposed to simulate the coupling to a 
complex reservoir with the help of Levy noise with vari- 
able exponent. In this manner, we extend usual quantum 
reservoir engineering in the spirit of the early studies of 
anomalous diffusion in complex media based on the ran- 
dom walk concept. The renewal process that we inves- 
tigate can actually be regarded as corresponding to a 
CTRW with power-law waiting time, but with a renor- 
malized diffusion range given by the localization length 
of the kicked rotator. 

At this point, it is appropriate to mention that anoma- 
lous diffusion has been shown before to occur in quantum 
kicked rotators subjected to deterministic aperiodic kick- 
ing. More precisely, subballistic wave packet spreading 
induced by quasiperiodic Fibonacci kick sequences has 
been established both away and at resonance [H, |40| . 
The effect of Levy noise at resonance has also been the 
subject of a recent study [4l|. An important difference to 
the present study is the fact that these studies do not con- 
sider the simultaneous presence of the periodic kicks in 
the intervals between the noise events. As a consequence, 
the classical momentum dynamics for the quasiperiodi- 
cally kicked systems is subdiffusive, and hence does not 
differ qualitatively from the quantum dynamics. The 
protocol proposed in the present paper stipulates the 
presence of such periodic kicks, which induce diffusive 



classical momentum dynamics, while quantum mechan- 
ically they drive the system towards localization. The 
ensuing qualitative difference between the noiseless quan- 
tum and classical dynamics indeed provides the terrain 
in which we explore the consequences of decoherence. 

The theory of decoherence that we have developed in 
this paper combines the properties of Levy noise, encoded 
in the generating function of the number of noisy events 
(jl4p . and a random-phase approximation, which exploits 
the complex quantum dynamics of the rotator appearing 
in the quasienergy transition amplitudes (|30p . It is worth 
emphasizing that it is nonperturbative (with respect to 
time) and applies both to stationary and nonstationary 
noise. 

One of the central results of our paper, contained in Eq. 
(|^T|) . is the demonstration of a regime of nonexponential 
loss of phase coherence when the noise is nonstationary 
(Levy exponent a < 1). This regime is characterized by 
an infinite coherence time, indicating that classical be- 
havior can never be attained, even though the amplitude 
of the noise is dramatically increased by several orders 
of magnitude. This remarkable property finds its origin 
in the divergence of the mean waiting time between ran- 
dom kicks and hence in the absence of a characteristic 
time scale in the engineered reservoir. 

The existence of a regime of slow decoherence has di- 
rect consequence for experimentally accessible physical 
quantities of the kicked rotator. We have first shown 
that the momentum spreading (j52p is subdiffusive and 
that the momentum profile retains a double-sided ex- 
ponential shape, see Fig. [31 These asymptotic features 
are in stark contrast to those in the usual decoherence 
regime, here realized for Levy exponent a > 1, which is 
characterized by normal diffusion and a Gaussian pro- 
file. The momentum distribution is routinely measured 
in kicked rotator experiments. In the atom-optical exper- 
iments @, 01 , the variance of the momentum is accessible 
up to times of the order of ten break times i*, indicat- 
ing that the subdiffusive dynamics depicted in Fig. [5] is 
measurable with current setups. 

Two other important quantities that we have consid- 
ered are the fidelity and the purity (or average fidelity), 
which are more traditional quantities for the character- 
ization of decoherence. Here again we have found clear 
signatures of the nonexponential decoherence regime. In 
Eq. ([M]) we have established a direct connection between 
the purity on the one hand and the decoherence factor 
and the inverse participation on the other. As a conse- 
quence, we could show that for a < 1 the decay of the pu- 
rity evolves from a stretched exponential to an algebraic 
decay, while a > 1 the decay changes from exponential to 
algebraic with exponent —1/2 (see Fig. [5|). The fidelity 
can be directly measured in echo experiments, see for 
example Ref. [36| , however to our knowledge such a mea- 
surement has not yet been performed in kicked rotator 
experiments. We also mention that we have found a di- 
rect relationship (|68p between the averaged logarithm of 
the fidelity and the sprinkling distribution of the renewal 
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process. 

Our concluding message is that reservoir engineering 
can be extended to mimic the coupling to complex non- 
stationary environments, allowing to control in a precise 
manner the loss of phase coherence and reveal hitherto 
unexplored decoherence scenarios. 
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